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Persistent currents in multichannel interacting systems 
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Persistent currents of disordered multichannel mesoscopic rings of spinless interacting fermions 
threaded by a magnetic flux are calculated using exact diagonalizations and self-consistent Hartree- 
Fock methods. The validity of the Hartree-Fock approximation is controled by a direct comparison 
with the exact results on small 4x4 clusters. For sufficiently large disorder (diffusive regime), the 
effect of repulsive interactions on the current distribution is to slightly decrease its width (mean 
square current) but to increase its mean value (mean current). This effect is stronger in the case 
of a long range repulsion. Our results suggest that the coupling between the chains is essential to 
understand the large currents observed experimentally. 

PACS numbers: 72.10.-d, 71.27.-|-a, 72.15.Rn 



The first observation of mesoscopic currents in very 
pure metallic nano-structures was done in pioneering ex- 
periments of Levy et al.^ and Chandrasekhar^'^. In the 
first case the experiment dealt with the average current of 
a system of 10^ disconnected rings in the diffusive regime 
while in the second a single ring was used. Although the 
existence of permanent currents in small metallic rings 
was predicted long ago'*'® the magnitude of the observed 
currents is still a real challenge to theorists. There has 
been a great activity to find a mechanism providing the 
correct magnitude of the currents. There is a general be- 
lief that the interaction has a crucial role in enhancing 
the current. But, so far, the role of the interaction in 
disordered systems is still unclear. Treating interaction 
and disorder on equal footings is a difficult task. 

In some previous work®'^ we have shown by exact di- 
agonalizations (ED) of small clusters that, for strictly 
ID systems of spinless fermions, the effect of a repulsive 
interaction is to increase further the localization of the 
electrons and hence to decrease the value of the current. 
Using a Hartree-Fock approach, Kato et al.*® have ob- 
tained qualitatively good agreement with our exact cal- 
culations. However, Giamarchi et al.® have shown that 
in the case of the ID Hubbard model (spin is included) 
the interaction enhances the persistent current. We then 
expect the increase of the current to be closely related to 
the decrease of the charge fiuctuations or, in other words, 
to the smoothening out of the charge density as this oc- 
curs in the ID Hubbard model whith repulsive interac- 
tion. Although this effect does not occur for ID spinless 
fermions, it might exist in higher dimensions. We shall 
therefore consider multi-channel rings relevant to exper- 
iments. Unfortunately in the case of multi-channel sys- 
tems the ED approach^ is limited to small rings (small 
number of channels) since the size of the Hilbert space 
increases exponentially fast with the number of lattice 
sites. In this paper, we use also an alternative and com- 
plementary approach based on a self-consistent Hartree- 
Fock (HE) treatment of the interaction between particles 



which enables one to treat much larger systems. Indeed, 
for increasing number of channels ie increasing dimen- 
sionnality quantum fiuctuations become less effective and 
a mean field treatment of the interaction might be appro- 
priate. However, we shall treat the quantum interference 
due to the disorder somehow exactly. Our method is 
different from the usual perturbative approach**^ where 
the corrections to the current due to the interacting term 
(Hartree and Eock terms) are calculated perturbatively. 
Recently, Ramin et al.*^ have numerically shown by such 
a perturbative approach that the first order correction 
to the second harmonic of the persistent current was 
in agreement with the analytical treatment**^; a near- 
est neighbour interaction tends to decrease the value of 
the typical current, while a repulsive extended Hubbard 
interaction enhances the typical current in any dimen- 
sion. On the contrary, in our self-consistent approach 
disorder and interaction are treated on equal footings. 
In another different approach Kopietz** by treating the 
Coulomb repulsion in a semi-classical way has found per- 
sistent currents in agreement with those measured by 
Chandrasekhar et al.^. 

This paper is organized as follow: first, the validity 
of the mean field treatment is controlled by comparing 
the HE results to the exact data obtained by Lanczos 
diagonalizations of small clusters. On the basis of this 
comparison an accurate estimation of the current in the 
HE approximation is suggested. Secondly, we use the 
HE method to handle larger clusters. We show that the 
interaction modifies the distribution of the currents. In 
the ballistic regime (small disorder) the interactions have 
little effects on the currents. However, in the diffusive 
regime. Coulomb repulsion between fermions can pro- 
duce a significant increase of the mean-current. We do 
not observe such a systematic effect in the case of a short 
range potential. 

The Hamiltonian reads: 
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where i, j are vectors of a 2D lattice (x,y) with Q < x < 
Lx — I and < y < Ly — I, Ly being the number of 
channels, Ly < L^- «i is the local density and Wi are on- 
site energies chosen randomly between -W/2 and W/2. 
In the following, we shall consider the cases of (i) a short 
range repulsion, Vij = V^i5(|i — j| — 1), and (ii) a Coulomb 
repulsion, Vij = — j|, i 7^ j. The hopping terms 

tij = i exp (2i7r$ij) are restricted to the bonds connect- 
ing nearest neighbor sites. Periodic boudary conditions 
are assumed in both directions. A flux $ (in units of 
$0 = hc/e) is applied through one hole of the 2D torus 
leading to a twist in the boundary conditions along x. 
This is realized e.g. by choosing $ij = -^Ux • (i — j) in 
(1). Hereafter, when not specified, energies are measured 
in unit of t. 

The current density operator is defined by = 
1 dn('S') 
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(exp (2i7r$/_Lj;) c| Ci+u, - h.c), 



(2) 



and the current in the ring is naturally given by /($) = 

ItF ( ^'ai^^ ) where Q means the expectation value in the 
ground state. Note that it can also be expressed alter- 
natively as /($) = ~ 27 ^ where E{<^) stands for the 
exact GS energy. 

Due to the random nature of the disorder, one can 
characterize the response to the fiux by the statistical 
distribution of the currents. Calculations are made for a 
fixed set of the parameters (V,W,(/)) and we average over 
at least 500 realisations of the disorder (w). The current 
distribution is investigated by calculating its first and 

second moments Ii = ^7($)^^ and I2 = \J (/^($))^ , 

where ( )^ means the average over disorder. Note that 
throughout the paper the moments of the currents are 
plotted in unit of 7o, average over half a period of 7^($) 
for V = Q and W = Q. The electron density is fixed to 
iri^ = 1/4 but we expect our results to be generic of the 
low density limit. 

As a first step, exact diagonalisations (ED) of (1) on 
a small 4x4 cluster are performed using a Lanczos 
algorithm^^ to calculate E(<^). Here we set $ = 1/4. 
Because of the limited length L^, = 4 of the ring we re- 
strict ourselves here to a repulsion between only nearest 
neighbor sites. In fig. 1(a) and (b) the first and second 
moments of the distribution of /($) are plotted versus 
V. We observe that the interactions do not enhance the 
current in this regime of disorder. Note that for such 
values of W the system is still in the ballistic regime. In- 
deed, the localization length is ^ ~ Lyl^, where Ly is the 
number of channels. If one estimates the elastic mean 
free path as /g ~ (2^7°)^ then for W/t < 4 one gets 
^ > Lx- A diffusive motion of the electrons is then diffi- 
cult to realize in such a small ring. On the contrary, for 



strong impurity scattering (localized regime) an increase 
of the currents was found^. These exact data will be- 
come very useful in the following as a necessary reference 
to control our approximation scheme. 

We now turn to the HF approximation of (1). The HF 
hamiltonian reads 
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which consists of a random potentiel (Hartree terms) 
w^'^'^ = Wi + Xlj^fi ^ij('^j) ^iid random hopping terms 
(Fock terms) tj/''^ = -|- Vij(cj Ci). The quantities (nj) 

and ^cj Ci^ are calculated self consistently. 

It is important to note that the HF current can be 
calculated alternatively by using /($) = ~27r'^^l$^ 
/($) = (Jj;). Although Ehf is the expectation value 
("W) in the HF groundstate which contains the fiux- 
dependance both in 7i and in the HF wavefunction, the 
self-consistency implies that (^^^jT^j^), where the vari- 
ation is made on the wavefunction only, is identically 
vanishing. A comparison between these two quantities 
provides then a useful check of the convergence of the 
self-consistent procedure. 

HF and exact results both obtained on the same 4x4 
cluster are compared in fig.l. We observe a rather good 
agreement between them and the same trends when V is 
applied. This provides then some evidence that the HF 
approach gives qualitatively correct results. 

We have observed for the 4x4 system a tendency of 
the localization length to decrease with increasing repul- 
sion in agreement with the fact that the current itself 
decreases. So far, our data suggest that the interaction 
do not enhance the current in the ballistic regime. This 
motivates the study of bigger systems where the diffusive 
regime can be more easily realized. 

^From now on we shall consider larger systems of size 
6x6 and 8x8 which can be handled only by using the 
HF approximation. It is important to increase simulta- 
neously Ly with Lx since, at fixed density, the number of 
energy level crossings increases with the number of chan- 
nels. We show in fig. 2 the current versus (j) for the clean 
system to illustrate this point. Level crossings lead to 
discontinuities of the current at some particular values of 
the fiux. The sign of the current depends on the parity 
of the number of electrons^® . We now compare the effect 
of short range and Coulomb repulsions in the presence 
of disorder. In fig. 3 and fig. 4 we ha ve plott ed Ii and 
the width of the distribution 61 = \/ Lj — If versus (j> . 
The calculations are done for W = 3, and V = and 0.5. 
Note that the sign of the current has not changed in com- 
parison with the clean case. First, we observe that the 
width of the distribution decrease with the interaction, 
the effect being larger in the case of short range repul- 
sion. Another important point is the large effect on the 
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first moment. We clearly observe that, in this diffusive 
regime, the Coulomb interaction enhances the current, 
while the screened interaction has a weaker effect. Note 
that the effect of the Coulomb repulsion is still small in 
the sense that it is far from cancelling out completely the 
drastic effect of the disorder potential to restore the clean 
current. However, we have to keep in mind that the rings 
considered here contain a small number of channels and, 
as shown in fig. 2, the number of level crossings is still 
quite small. Considering this last point we argue that 
the effect observed here is in fact quite significant. We 
expect the effect to become much larger when the num- 
ber of channels (and hence of level crossings) becomes 
larger. Note that similar increase of the currents have 
also been obtained recently by exact diagonalizations of 
4x4 clusters in the large W regime^ 

In conclusion, the HF results are in qualitatively good 
agreement with the exact calculations. We find that, in 
the case of spinless fermions, the Coulomb interaction 
enhances the current, while a screened interaction has a 
weaker effect. The increase (decrease) of the current is al- 
ways consistently related to the increase (decrease) of the 
localization length i.e. to a delocalization (localization) 
of the electronic wavefunctions in average. Qualitatively 
we expect a fraction of the Fermi gas to become local- 
ized at the minima of the potential and, because of the 
interaction, to screen out these minima for the rest of 
the electrons. This physical picture is particularly clear 
in the HF equations themselves where a new effective 
potential is generated by the interaction. Whether suffi- 
cient screening can be self-consistently generated depends 
on the nature of the interaction between fermions. Our 
work strongly suggests that a long range interaction is 
more effective to smooth out the impurity potential and 
hence to increase the current. We also argue that this 
effect will become stronger and stronger as the number 
of channels is increased. A detailed study of the role of 
the transverse dimension is left for future work. 
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FIG. 1. First moment /i (full symbols) and second moment 
I2 (open symbols) of the current distribution at 1/1 = 1/4 for 
W=l (a) and W=3 (b). Squares (triangles) correspond to the 
HF (exact) current ( Jx). The average is taken over Nr = 500 
realisations of the disorder in the ED case and over Nr = 1000 
in the HF approximation. Error bars of order 81 1 Nr are in- 
cluded unless smaller than the size of the dots. 

FIG. 2. First moment /i versus 4> aX V = Q and IF = for 
4x4, 6x6 and 8x8 systems. 

FIG. 3. Moments of the current of a 6x6 site system for 
IF = 3 . Each symbol represents an average over 500 con- 
figurations of the disorder, (a) First moment /i for = 
(triangles), a short range interaction = 0.5 (squares) and 
a Coulomb repulsion V = 0.5 (pentagons), (b) width of the 
distribution 81 for = (triangles), a short range interac- 
tion y = 0.5 (squares) and a Coulomb interaction = 0.5 
(pentagons). Typically, the error bars corresponding to the 
I\ data (not shown) are of the order of ^/2. 



FIG. 4. Same as fig. 3 for a 8x8 site system. 
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